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ABSTRACT 

The  PSR  B1259-63  system  (Johnston  et  al.  1992,  1994)  was  observed  near 
periastron  by  the  Compton  Gamma-Ray  Observatory  in  January  1994.  This 
system  contains  a  rapidly  rotating  pulsar  and  a  Be  star  in  a  highly  eccentric 
binary  orbit.  We  report  the  discovery  by  the  OSSE  instrument  of  unpulsed 
emission  with  a  hard  power-law  spectrum  between  50  and  200  keV  from  the 
direction  of  this  system.  Neither  diffuse  galactic  background  emission  nor 
nearby  X-ray  binaries  contribute  significantly  to  the  detected  flux.  Our 
results  are  particularly  important  for  the  theory  of  interaction  of  pulsars  with 
gaseous  environments.  We  interpret  the  hard  X-ray  emission  as  synchrotron 
radiation  from  relativistic  particles  of  the  PSR  B1259-63  wind  being  shocked  and 
accelerated  within  the  binary.  Our  results  indicate,  for  the  first  time  in  a  binary 
pulsar,  that  shock  acceleration  can  increase  the  original  energy  of  pulsar  wind 
particles  by  a  factor  ~  10,  despite  the  high  synchrotron  and  inverse  Compton 
cooling  rates  near  periastron.  The  derived  shock  properties  (efficiency,  radiation 
spectrum,  timescale)  are  relevant  for  a  broad  class  of  high-energy  astrophysical 
sources  characterized  by  shocked  relativistic  plasmas  subject  to  strong  radiative 
cooling. 
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Subject  headings:  gamma  rays:  observation  —  stars:  pulsars:  individual 
(PSR  B1259-63)  —  shock  waves  —  acceleration  of  particles 


1.  Introduction 


The  PSR  B1259-63  system  contains  a  rapidly  rotating  radio  pulsar  with  spin  period 
P  =  47.76  ms  and  spindown  luminosity  —  9  x  1035  erg  s  1  orbiting  around  a  massive 
Be  star  companion  in  a  highly  elliptical  orbit  with  orbital  period  ~  3.4  yrs  (Johnston 
et  al.  1992,  1994).  The  interest  in  a  detailed  study  of  high-energy  emission  from 
the  PSR  B1259-63  system  near  periastron  is  manifold.  Be  stars  are  known  to  produce 
substantial  gaseous  outflows  mainly  concentrated  in  their  equatorial  planes  (e.g.,  Waters 
et  al.  1988),  and  the  interaction  of  nebular  material  with  a  radio  pulsar  can  in  principle 
lead  to  accretion-powered  high-energy  emission  or  detectable  shock  emission.  High-energy 
emission  from  the  PSR  B1259-63  system  is  likely  to  be  enhanced  near  periastron,  since  the 
pulsar  approaches  its  Be  star  companion  at  an  orbital  distance  a  ~  1013  cm,  corresponding 
to  ~  10  —  20  stellar  radii  of  the  companion.  Typical  orbital  distances  near  apastron  of 
the  PSR  B1259-63  system  are  a  factor  of  ten  larger.  Radio  pulsars  lose  their  rotational 
energy  by  the  emission  of  an  electromagnetic  and  relativistic  particle  wind  composed  of 
electrons  and  positrons  (Kennel  &  Coroniti  1984).  It  is  also  likely  that  the  relativistic 
winds  from  radio  pulsars  contain  heavy  ions  (see  Gallant  &  Arons  (1994)  and  references 
therein  for  evidence  supporting  this  view).  In  the  PSR  B1259-63  system  the  interaction 
of  a  relativistic  pulsar  wind  and  a  mass  outflow  from  a  Be  star  can  be  studied  for  the  first 
time  in  a  time- variable  environment  with  a  shock  radius  rs  ~  1013  cm  (Tavani,  Arons,  & 
Kaspi  1994). 

High-energy  emission  (i.e.  X-rays  and  gamma-rays)  can  be  radiated  by  a  pulsar  wind 
termination  shock  during  the  periastron  passage  of  PSR  B1259-63,  and  this  possibility 
motivated  a  long  (three-week)  observation  by  the  Compton  Gamma-Ray  Observatory 
(CGRO).  The  CGRO  instruments  cover  a  broad  energy  range  between  30  keV  and  ~30  GeV. 
Here  we  focus  on  the  positive  detection  of  hard  X-rays  in  the  direction  of  the  PSR  B1259-63 
system  by  the  OSSE  instrument.  We  briefly  comment  on  how  the  detected  emission  can  be 
used  to  sensibly  constrain  models  of  high  energy  emission  from  pulsar  termination  shocks. 
A  more  detailed  theoretical  analysis  of  our  results  will  appear  elsewhere  (Tavani  &  Arons 
1995). 


2.  Observations 


OSSE  observed  PSR  B1259-63  for  the  three-week  period  3-23  January  1994  chosen  to 
cover  the  periastron  passage  (9  January  1994,  Johnston  et  al.  1994).  The  OSSE  instrument 
consists  of  four  nearly  identical  large-area  NaI(Tl)-CsI(Na)  phoswich  detector  systems 
(Johnson  et  al.  1993).  It  covers  the  energy  range  from  50  keV  to  10  MeV  with  good 
spectral  resolution.  A  tungsten  slat  collimator  defines  the  3.8°  X  11.4°  field  of  view,  which 
was  chosen  as  a  compromise  allowing  high  sensitivity  to  both  diffuse  emission  and  point 
sources.  Spectra  are  accumulated  in  a  sequence  of  two-minute  measurements  of  the  source 
field  alternated  with  two-minute,  offset-pointed  measurements  of  background. 

The  pulsar  is  located  in  the  galactic  plane  at  304.2°  longitude,  where  galactic  diffuse 
continuum  emission  is  detectable  by  OSSE.  To  minimize  the  effect  of  the  galactic  emission, 
the  observation  of  PSR  B1259-63  was  performed  with  the  long  axis  of  the  collimator 
oriented  perpendicular  to  the  plane  and  centered  on  the  pulsar.  Background  fields  on  either 
side,  centered  at  308.7°  and  296.4°  longitude,  were  selected  to  avoid  nearby  point  sources 
such  as  the  X-ray  binaries  GX301-2,  Cen  X-3,  and  2S  1417-624.  These  objects,  shown  by 
the  BATSE  instrument  on  CGRO  to  be  occasional  emitters  in  the  20-50  keV  band  (R. 
Wilson,  private  communication),  are  each  suppressed  by  a  factor  of  >  102  by  the  collimator 
and  account  for  no  detectable  emission  in  the  source  field.  Unavoidably,  the  source  field 
contained  the  Be-star  binary  pulsar  GX304-1;  however,  this  X-ray  pulsar  has  not  been 
detected  by  BATSE  in  three  years  of  daily  observation  (M.  Finger,  private  communication), 
and  it  was  shown  by  EXOSAT  to  be  in  an  off  state  in  1984  (Pietsch  et  al.  1986).  A  total 
exposure  time  T  =  1.9  X  105  s  of  highest-quality  data  were  collected  in  two  OSSE  detectors 
on  the  PSR  B1259-63  source  field,  with  approximately  an  equal  time  on  the  background 
fields. 


3.  Results 


The  energy  spectrum  from  the  PSR  B1259-63  region  is  shown  in  Fig.  1.  Emission  is 
detected  from  50  keV  to  ~200  keV,  with  a  total  statistical  significance  of  4.8cr,  at  a  level 
of  ~5  mCrab  flux  units. 

The  best-fit  power  law  photon  number  spectrum  for  photon  energies  E  in  the  50-300 
keV  band  is  dN1/dE  =  (2.8  ±  0.7)  x  10 7  cm-2  s-1  MeV-1,  with  photon  index 
a  =  1.8  ±  0.6  and  E2  =  E'/(100  keV).  The  power-law  fit  is  good  (%|of  =  0-90  f°r  40  degrees 
of  freedom).  Also  shown  in  Fig.  1  are  schematic  extrapolations  from  the  power-law  fits  to 
the  ASCA  data  from  28  December  1993,  10  January  1994,  and  26  January  1994  (Kaspi 
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et  al.  1995).  We  note  that  the  uncertainties  in  ASCA’s  fit  parameters  are  such  that  the 
three  models  are  essentially  equivalent  at  100  keV.  The  emission  in  the  OSSE  band  is 
therefore  consistent  in  intensity  and  spectral  shape  with  a  simple  power-law  extrapolation 
of  the  ASCA  observations.  The  continuation  of  the  power  law  spectrum  to  hard  X-ray 
energies  is  very  significant  for  the  interpretation  of  the  emission  (see  Discussion).  At  2.0  kpc 
(Johnston  et  al.  1994),  the  inferred  50-200  keV  luminosity  of  the  PSR  B1259-63  system  is 
Lx  —  3  X  1034  erg  s-1,  and  the  hardness  of  the  OSSE  spectrum  suggests  that  the  luminosity 
per  energy  decade  peaks  near  or  above  200  keV. 

There  is  no  evidence  for  variability  in  the  OSSE  data;  however,  because  of  the  relatively 
low  statistical  significance  of  the  detection,  the  constraints  on  variability  are  not  stringent. 
The  95%-confidence  upper  limits  on  50-200  keV  variability  on  daily  or  weekly  timescales  are 
~190%  and  ~70%,  respectively.  Any  flux  decrease  by  ~50%  during  the  3-week  observing 
period,  as  suggested  by  the  ASCA  data,  is  undetectable  by  OSSE.  Because  of  the  brevity 
of  our  observation  in  relation  to  the  3.4-year  binary  period,  we  are  unable  to  place  a 
meaningful  limit  on  modulation  at  orbital  timescales. 

For  weak  detections  such  as  that  from  the  PSR  B 1259-63  field,  it  is  important  to 
address  possible  systematic  errors  or  alternate  sources  for  the  emission.  A  study  of 
residuals  in  the  50-200  keV  band  from  an  aggregrate  of  many  OSSE  pointings  indicates 
that  the  systematic  errors  are  smaller  than  the  statistical  errors  and,  therefore,  that  the 
~5cr  excess  reported  here  should  be  exceedingly  rare.  The  background  fields  of  view  were 
chosen  to  exclude  known  point  sources  of  hard  X-rays  in  the  region,  and  the  source  field 
otherwise  contained  only  GX304-1,  which  is  apparently  in  an  off  state  (see  Observations). 
Furthermore,  we  find  no  evidence  for  pulsations  from  GX304-1  in  the  OSSE  data  at  periods 
between  267  s  and  275  s,  which  reasonably  bound  extrapolations  of  the  historical  period. 

OSSE  is  conducting  an  extensive  survey  of  the  galactic  diffuse  line  and  continuum 
emission  (e.g.  Purcell  et  al.  1994).  To  date,  measurements  appropriate  to  estimating 
the  diffuse  emission  in  the  PSR  B1259-63  region,  i.e.  those  at  similar  separation  from 
the  galactic  center  and  without  significant  contributions  from  known  point  sources,  are 
available  only  on  the  opposite  side  of  the  galactic  center,  at  longitudes  40°  and  58°.  From 
these  measurements  and  a  simple  one- dimensional  galactic  plane  model,  we  estimate  that 
the  large-scale  (i.e.  ~  10°)  longitudinal  gradient  in  the  emission  per  degree  in  the  50-200 
keV  band  is  ~  —3  X  10-5  7  cm-2  s-1  MeV-1  deg-2.  Assuming  that  the  diffuse  emission 
is  generally  symmetric  about  the  galactic  center  and  convolving  this  emission  through  the 
instrument  response,  we  estimate  that  the  residual  diffuse  50-200  keV  flux  for  our  viewing 
strategy  of  the  PSR  B1259-63  region  is  ~  2  X  10-4  7  cm-2  s-1  MeV-1,  a  factor  of  ~10  less 
than  the  observed  flux.  We  note,  however,  that  the  level  of  small-scale,  local  fluctuations  in 
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the  diffuse  emission  is  somewhat  uncertain  and  not  easily  measured,  given  the  presence  of 
nearby  point  sources.  Such  fluctuations  are  the  most  likely  alternative  explanation  for  the 
observed  flux. 

Using  the  radio  ephemeris  (Manchester  et  al.  1995)  for  the  spin  and  binary  orbital 
parameters  of  PSR  B1259-63  ,  we  searched  for  pulsed  emission  at  the  rotation  period  of  the 
neutron  star  in  OSSE’s  high  time  resolution  data  channels,  which  provided  8-ms  (or  faster) 
rate  samples  over  50-180  keV.  No  pulsed  emission  was  detected  on  any  day  or  for  the  total 
observation.  However,  we  caution  that  because  of  the  long  integrations  required  in  the  hard 
X-ray  regime,  small  errors  in  the  pulsar  ephemeris  could  smear  coherent  pulsations.  The 
95%-confidence  upper  limits  for  pulses  of  phase-width  0.15  or  0.5  in  the  50-180  keV  band 
over  the  21-day  observation  period  are  $0.is  =  3  X  10-4  and  $0.s  =  7  X  10-4  7  cm-2  s-1 
MeV-1,  respectively.  The  former  corresponds  to  3  X  10-3  Crab  pulsar  flux  units.  Upper 
limits  for  individual  days  are  typically  a  factor  of  \/21  higher. 

The  region  of  PSR  B1259-63  was  also  observed  by  OSSE  in  1992  September  19  -  1992 
October  3  (TJD  8518  -  8532),  when  the  pulsar  was  far  from  periastron,  as  part  of  OSSE’s 
galactic  plane  survey.  However,  the  sensitivity  was  reduced  by  the  placement  of  the  pulsar 
in  the  wing  of  the  collimator,  which  was  aligned  with  the  plane.  Any  unpulsed  emission 
from  the  pulsar  was  masked  by  GX301-2  and  Cen  X-3  (both  at  high  sensitivity  near  the 
center  of  the  field  of  view)  as  well  as  by  diffuse  emission  from  the  plane.  While  these  data 
therefore  cannot  provide  a  useful  limit  on  unpulsed  emission  from  PSR  B1259-63  away  from 
periastron,  they  have  been  used  in  a  search  for  pulsed  emission  (Ray  et  al.  1992).  The 
resulting  99.9%-confidence  limit  of  6  X  10-3  7  cm-2  s-1  MeV-1  in  the  64-150  keV  band  was 
achieved  through  a  search  of  a  large  volume  of  acceleration  space,  as  an  adequate  binary 
ephemeris  was  not  then  available  in  the  radio.  We  have  reanalyzed  this  data  with  the 
ephemeris  of  Johnston  et  al.  (1994)  and  set  an  upper  limit  of  1.4  X  10-3  7  cm-2  s-1  MeV-1. 


4.  Discussion 


Hard  X-ray  emission  from  a  neutron  star  with  a  power-law  extending  up  to  ~  200  keV 
and  with  photon  index  a  ~  1.8  —  2  is  a  manifestation  of  energetic  particle  acceleration. 
Accreting  neutron  stars  do  not  usually  emit  in  the  hard  X-ray  range  (White  et  al.  1995) 
with  photon  energies  Ex  ~  50  keV.  Only  sporadic  hard  X-ray  emission  from  a  few  accreting 
neutron  stars  has  been  reported  (Barret  &  Vedrenne  1994),  and  a  power-law  fit  to  the 
spectrum  typically  gives  a  photon  index  ~  2  —  3  or  larger.  In  the  case  of  emission  from  the 
direction  of  the  PSR  B1259-63  system,  the  lack  of  pulsations,  the  absence  of  detectable 
flux  variations  within  the  timescale  of  weeks,  the  relatively  low  hard  X-ray  luminosity,  and 
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the  spectrum  substantially  harder  than  for  accreting  neutron  stars  (and  consistent  with  an 
extrapolation  to  soft  X-ray  energies)  all  argue  against  accretion-powered  emission  (Tavani 
&  Arons  1995,  Kaspi  et  al.  1995). 

Hard  X-rays  provide  crucial  information  on  the  shock  emission  due  to  the  interaction 
of  the  relativistic  particle  wind  of  PSR  B1259-63  with  the  nebular  gas  of  the  Be  star  mass 
outflow.  As  in  the  case  of  the  Crab  nebula  (Kennel  &  Coroniti  1984,  Gallant  &  Arons  1994), 
shock  acceleration  and  synchrotron  radiation  capable  of  producing  unpulsed  hard  X-ray 
emission  occurs  at  a  shock  radius  within  the  PSR  B1259-63  binary  where  pressure  balance 
is  established  (Tavani  et  al.  1994)  between  the  pulsar  relativistic  wind  (electrons/positrons, 
possibly  ions)  and  the  ram  pressure  of  the  Be  star  outflow.  The  radiative  environment  near 
the  shock  radius  within  the  PSR  B1259-63  system  is  quite  different  from  that  of  the  Crab 
nebula,  and  the  large  synchrotron  and  inverse  Compton  cooling  rates  (~  10-2  —  10-3s-1) 
at  the  PSR  B1259-63  periastron  make  the  detection  of  high-energy  emission  a  unique 
diagnostic  of  shock  acceleration  subject  to  strong  radiative  cooling  (Tavani  et  al.  1994, 
Tavani  &  Arons  1995).  By  comparing  the  OSSE  and  ASCA  results  a  unified  picture  for  the 
emission  and  spectrum  emerges.  Spectral  calculations  for  models  with  a  weak  ram  pressure 
of  the  Be  star  outflow  show  that  the  dominating  inverse  Compton  spectrum  is  typically 
much  too  flat  to  agree  with  the  OSSE  data  (Tavani  &  Arons,  1995).  Instead,  the  spectral 
shape  and  intensity  of  the  power-law  emission  from  the  PSR  B1259-63  system  are  consistent 
with  synchrotron  radiation  of  shock-accelerated  electron-positron  pairs  for  an  intermediate 
value  of  the  ram  pressure  of  the  Be  star  outflow  and  a  shock  distance  relatively  close  to  the 
pulsar  rather  than  to  the  Be  star.  The  inferred  characteristics  of  the  electron-positron  pairs 
in  the  pulsar  wind  are  a  Lorentz  factor  of  71  ~  106  and  a  ratio  of  electromagnetic  energy 
density  to  particle  kinetic  energy  density  in  the  wind  of  a  ~  10-1  —  10-2  (Tavani  &  Arons, 
1995).  The  efficiency  of  conversion  of  the  pulsar  wind  kinetic  and  electromagnetic  energy 
into  radiation  in  the  50-200  keV  band  is  ~  3%  for  a  system  distance  of  2  kpc  (Johnston  et 
al.  1994).  The  absence  of  a  spectral  break — from  the  soft  X-ray  band  above  em8„  ~  1  keV 
(Kaspi  et  al.  1995)  through  the  lower  limit  on  the  high-energy  cutoff,  emax  ~  200  keV,  from 
the  OSSE  observations — strongly  constrains  the  acceleration  mechanism.  The  collisionless 
pulsar  wind  shock  must  produce  a  power-law  particle  distribution  function  of  index  5 
between  a  minimum  e±-pair  energy  E\  =  71  me  c2  and  an  upper  energy  cutoff  Ec  =  7cme  c2, 
with  me  the  electron-positron  mass  and  c  the  speed  of  light.  For  a  synchrotron  model  of 
shock  emission,  we  obtain 

„  10  5  (i) 

'Jl  x  ^ min  / 

and  6  =  2a  —  1  ~  2  —  3.  The  OSSE  observation  of  the  PSR  B1259-63  system  therefore 
demonstrates,  for  the  first  time  in  a  galactic  source,  the  ability  of  shock  acceleration  to 
efficiently  energize  a  relativistic  plasma  subject  to  radiative  cooling  with  typical  timescale 
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tc  ~  102  s.  The  important  derived  constraints  on  the  radiation  efficiency  (~  3%), 
acceleration  timescale  rc),  index  of  the  post-shock  particle  distribution  (6  ~  2  —  3)  are 
of  crucial  interest  for  any  theory  of  relativistic  shock  acceleration. 

Our  results  can  have  a  broad  range  of  applications  to  astrophysical  sources  with 
relativistic  plasmas  being  shocked  and  accelerated  in  a  strongly  radiative  environment. 
Our  observations  of  the  PSR  B1259-63  system  show  that  shock  emission  is  an  efficient 
mechanism  of  unpulsed  high-energy  radiation  in  pulsar  binaries.  Poorly  understood 
unidentified  and  time- variable  high-energy  galactic  sources  (e.g.,  Fichtel  et  al.  1994)  may 
reveal  energetic  pulsars  powering  unpulsed  emission  in  binaries. 
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Fig.  1. —  OSSE  spectrum  of  emission  from  the  region  of  PSR  B1259-63.  Schematic 
extrapolations  to  the  power-law  spectra  of  ASCA  observations  from  1993  December  28  (solid 
line),  1994  January  10  (dotted  line),  and  1994  January  26  (dashed  line)  are  also  shown.  ASCA 
extrapolations  obtained  from  the  analysis  of  Kaspi  et  al.  (1995). 


